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ABSTRACT 

We compute the emission of gravitational radiation from the merging of a close white 
dwarf binary system. This is done for a wide range of masses and compositions of 
the white dwarfs, ranging from mergers involving two He white dwarfs, mergers in 
which two CO white dwarfs coalesce to mergers in which a massive ONe white dwarf 
is involved. In doing so we follow the evolution of binary system using a Smoothed 
Particle Hydrodynamics code. Even though the coalescence process of the white dwarfs 
involves considerable masses, moving at relatively high velocities with a high degree 
of asymmetry we find that the signature of the merger is not very strong. In fact, the 
most prominent feature of the coalescence is that in a relatively small time scale (of 
the order of the period of the last stable orbit, typically a few minutes) the sources 
stop emitting gravitational waves. We also discuss the possible implications of our 
calculations for the detection of the coalescence within the framework of future space- 
borne interferometers like LISA. 
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1 INTRODUCTION 

Gravitational waves are a direct consequence of the General 
Theory of Relativity. Many efforts have been done so far to 
detect them, but due to the intrinsic experimental difficulties 
involved in the process of detection and in the data analysis 
no definite result has been yet obtained. Supernova core col- 
lapse, binary systems involving compact objects, like double 
black holes, double neutron stars or double white dwarfs are, 
amongst others, promising sources of gravitational waves — 
see Schutz (1999) for a comprehensive review on the subject. 
Moreover, with the advent of the current generation of ter- 
restrial gravitational wave detectors, like LIGO (Abramovici 
et al. 1992), VIRGO (Acernese et al. 2004), GEO600 (Willke 
et al. 2004), or TAMA (Takahashi et al. 2004), and of space- 
borne interferometers like LISA 1 (Bender et al. 1998, 2000), 
gravitational wave astronomy will probably be soon a tan- 
gible reality. 

As already mentioned, one of the most promising 
sources of gravitational waves are galactic binary systems 
containing at least one compact object. Galactic bina- 
ries, such as neutron stars binaries, cataclysmic binaries or 
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close white dwarf binaries, are guaranteed sources for LISA 
(Mironowsky 1965; Evans, Iben & Smarr 1987), provided 
that the sources are at sufficiently close distances. In fact, 
emission from galactic close white dwarf binary systems is 
expected to be the dominant contribution to the background 
noise in the low frequency region, which ranges from ~ 10 -3 
up to ~ 10~ 2 Hz (Bender et al. 1998). Additionally, from 
very simple considerations about the initial mass funtion, it 
is easy to see that galactic close white dwarf binaries must 
be quite common (Hils, Bender & Webbink 1990) and, con- 
sequently, if the amplitude of the gravitational waves is large 
enough we should be able to eventually detect them during 
the operation of LISA. Moreover, the merging of two white 
dwarfs by emission of gravitational radiation will be the fi- 
nal destiny of a good fraction of this type of binary systems. 
Since during the merging process a sizeable amount of grav- 
itational waves is expected to be produced (Guerrero et al. 
2004) it is important to characterize which would be the 
gravitational wave emission of such process and to assess 
the feasibility of dectecting them. 

The process of formation of close white dwarf bina- 
ries involves two mass transfer episodes of the progenitor 
stars when each of the components of the binary system 
evolves off the main sequence. Depending on when dur- 
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ing the lives of the binary components the mass transfer 
episodes occur the components may have different core com- 
positions. In particular, we may have He-He systems with a 
total mass Mtot < 0.75 Mq, He-CO for those systems with 
masss within the range 0.75 Mq < M tot < 1.45Mq, CO-CO 
for masses larger than Mtot ~ 1-45 Mq and even He-ONe or 
CO-ONe systems when one of the white dwarfs is a massive 
one. Although the astrophysical scenarios in which a merger 
of two white dwarfs in a close binary system can occur and 
their relative frequencies have been relatively well studied 
— see, for instance, Yungelson et al. (1994), and Nelemans 
et al. (2001a, 2001b), and references therein — the process 
of merging itself has received little attention until very re- 
cently. Indeed, one of the probable reasons for this lack of 
theoretical models is the heavy computational demand in- 
volved in the simulation of an intrinsically three-dimensional 
phenomenon. However, in sharp contrast, the coalescence of 
two neutron stars has been extensively studied — see, for 
instance, Rosswog et al. (2000), Rosswog & Davies (2002), 
and Rosswog & Liebendorfer (2003), and references therein, 
for some of the most recent works on this subject. 

In a recent paper (Guerrero et al. 2004) we thor- 
oughly studied the merging of white dwarf binary sys- 
tems for a wide range of masses and compositions. In do- 
ing this we used an up-to-date Smoothed Particle Hydro- 
dynamics code. This method was first proposed by Lucy 
(1977) and, independently, by Gingold & Monaghan (1977). 
The fact that the method is totally Lagrangian and does 
not require a grid makes it specially suitable for study- 
ing an intrinsically three-dimensional problem like the co- 
alescence of two white dwarfs. In the present paper we 
will not discuss the simulations presented there. The inter- 
ested reader can find them in mpeg format at the follow- 
ing URL: ftp://ftp.ieec.fcr.es/pub/astrofisica/SPH 
Instead, in this paper we compute the gravitational wave 
signature expected from the merging of two white dwarfs as 
obtained from those SPH calculations. The paper is orga- 
nized as follows. In §2 we briefly describe the main ingre- 
dients of our SPH code, whereas in §3 we summarize the 
physics of the emission of gravitational waves. In §4 we dis- 
cuss the calibration and consistency checks we have done in 
order to assess the reliability of our results, whereas in §5 
and §6 we present our results. Finally, in §7 we summarize 
our results and draw our conclusions. 



2 THE SMOOTHED PARTICLE 
HYDRODYNAMICS CODE 

We follow the coalescence of the binary system using a La- 
grangian particle numerical code and, more specifically, a 
Smoothed Particle Hydrodynamics (SPH) code. Our SPH 
code has been described at length in Guerrero et al. (2004) . 
However, for the sake of completeness we provide here a 
short description of the main input physics and numeri- 
cal techniques. Our code follows closely the prescriptions 
of Benz (1990), where the basic numerical scheme for solv- 
ing the hydrodynamic equations can be found. We use 
the standard polynomic kernel of Monaghan & Lattanzio 
(1985). The gravitational forces are evaluated using an oc- 
tree (Barnes & Hut 1986). The artificial viscosity adopted in 
this work is that of Balsara (1995). Regarding the integra- 



tion method we use a predictor-corrector numerical scheme 
with variable time step (Serna, Alimi & Chieze 1996) , which 
turns out to be quite accurate. The adopted equation of state 
is the sum of three components. The ions are treated as an 
ideal gas but taking into account the coulombian corrections. 
We have also incorporated the pressure of photons. Finally, 
the most important contribution is the pressure of degen- 
erate electrons which is treated as the standard zero tem- 
perature expression plus the temperature corrections. The 
nuclear network adopted here (Benz, Hill, & Thielemann 
1989) incorporates 14 nuclei: He, C, O, Ne, Mg, Si, S, Ar, 
Ca, Ti, Cr, Fe, Ni and Zn. The reactions considered are cap- 
tures of a particles, and the associated back reactions, the 
fusion of two C nuclei, and the reaction between C and O 
nuclei. All the rates are taken from Rauscher & Thielemann 
(2000). The thermal evolution of the system is followed in 
two ways. On the one hand the variation of the internal 
energy is followed according to: 
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where W(rij,h) is the smoothing kernel and the rest of 
the symbols have their usual meaning. The smooting length 
adopted in in this work is the linear average of the par- 
ticles considered in the corresponding calculation; that is 
hij = (hi + hj) /2,On the other, we simultaneously compute 
the temperature variation according to: 
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If in the region under study the temperature is below 6 x 
10 8 K or the density is smaller than 6 x 10 3 g/cm 3 we use 
equation (2) otherwise we use equation (1). We have found 
that in this way energy is better conserved. 



3 EMISSION OF GRAVITATIONAL WAVES 

We compute the gravitational wave emission in the slow- 
motion, weak-field quadrupole approximation (Misner et al. 
1973). The dimensionless wave strain, h, in the transverse- 
traceless gauge is given by: 



<_tt / x ^ 2Gd 2 Q]?(t-R) 



h jk (t, x) = 



c 4 d 
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where t — R = t — d/cis the retarded time, d is the distance 
to the observer and Qj^(t — R) is the reduced quadrupole 
moment of the mass distribution, which is given by: 

Qjk(t - R) = J p(x, t - R)(x 3 x k - ^x 2 S jk )d 3 x (4) 

The rest of the symbols have their usual meaning. It is useful 
to express the time derivative of the quadrupole moment in 
the following way (Nakamura & Oohara 1989): 

Qjk(t-R) = Pi jk i(N) J d 3 xp[2v k v l -x k a l <j>- x l d k cj>j (5) 
where 
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is the transverse-traceless projection operator onto the plane 
ortogonal to the outgoing wave direction, N, and <f> is the 
gravitational potential. Now, one can express Eq. © in the 
following way: 



(7) 



hjk (t,x) = -j^( Jl 4+(t,x)e +3 -fc + J 4 x (t,x)e X3 i 



where the polarization tensor coordinate matrices are de- 
fined as: 
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the dimensionless amplitudes h+ = A+/d and h x = A x /d 
are the two independent modes of polarization in the 
transverse-traceless gauge, and the amplitudes are respec- 
tively given by 

A+(t, x) = Q xx - Q yy , A x (t, x) = +2Q xy (9) 
for i = 0, and 

A+(t,x) =Q ZZ -Q yv , A x (t,x) = -2Q yz (10) 
for i = 7r/2. 

In our case we have a collection of n individual SPH 
particles. Consequently, Eq. must be discretized and it 
is computed according to the following expression: 



QJ^t-R) « P ljfc; (N)^m(p)[2v fc (p)v i (p) 

p=l 

+ x*(p)a , (p)+x , (p)a*(p)] (11) 

Where m(p) is the mass of each SPH particle, and x(p), 
v(p) and a(p) are, respectively, its position, velocity and 
acceleration. 



4 CALIBRATION AND CONSISTENCY 
CHECKS 

We have done two tests, gravitational wave emission from 
a single, isolated star, and gravitational wave emission from 
a close white dwarf binary system in a circular orbit. For 
both cases there exist analytical solutions to which we can 
compare our numerical results. In the first case, we have 
followed the time evolution of an isolated 1 Mq star using 
2 x 10 4 SPH particles of the same mass. For the second test 
we have followed the evolution of a binary system made of 
two white dwarfs of the same mass (1 Mq) in a circular 
orbit. Each one of the white dwarfs was simulated using 
2 x 10 4 SPH particles of the same mass. 

The first of our tests was designed to set the zero point 
of our calculations. Since in SPH simulations the particles 
are allowed to move freely under the action of their own 
gravitational potential and of the pressure forces, and since 
the mass of each particle is relatively large it is not obvi- 
ous a priori whether or not a stable relaxed configuration 
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Figure 1. Gravitational wave emission from an isolated, spherical 
star. The dimensionless strains and h x are measured in units 
of 10~ 25 . The source is located at a distance of 10 kpc. 



radiates gravitational waves. Figurc^shows the dimension- 
less strains for the case of an isolated white dwarf. We only 
show times larger than 60 s, for which the star is already 
relaxed to its final configuration. The relaxation procedure 
consisted in allowing the initial configuration (which con- 
sisted in randomly distributing the SPH particles according 
to the density profile of a zero temperature white dwarf of 
the same mass) to evolve for a long enough time until the os- 
cillations of the resulting configuration were completely neg- 
ligible. In this way we check whether or not the numerical 
noise produces a negligible emission of gravitational waves. 
And this is indeed the case. As it can be seen the emission 
of gravitational waves is negligible, as it should be, given 
that the relaxed configuration presents spherical symmetry 
— see Eq. 

With regard to the second test, it is worth recalling that 
the emission of gravitational waves from a binary system, 
can be obtained quite easily by assuming that both stars 
are point-like mass distributions. By doing so, from Eqs. ^ 
and one obtains 



h. 



(1 + cos i) cos 2ut 



and 



h > 



2V2^^-(wM tot )t(cosi) sin 2wt, 



(12) 



(13) 



where uj is the angular velocity of the stars, is the re- 
duced mass, M is the mass of the white dwarfs, and i is the 
observation angle with respect to the orbital plane. 

In our SPH simulations we have chosen the binary 
system to have a separation of 0.05 Rq or, equivalently, 
uj — 7.94 x 10 -2 s _1 . The two white dwarfs describe circular 
orbits and no mass is transferred between both components. 
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Figure 2. Gravitational wave emission from a close white dwarf 
binary system. The dimensionless strains and h x are mea- 
sured in units of 10 — 22 . The solid line corresponds to the theo- 
retical solution whereas the dots correspond to our numerical so- 
lution. Again, the source is assumed to be at a distance of 10 kpc. 



Figure 3. Residuals between the theoretical solution and the 
gravitational wave emission from the close white dwarf binary 
system of Fig. [5] Note that the dimensionless strains h+ and h x 
are measured in units of 10 — 23 , one order of magnitude smaller 
than the scale of Fig. |21 



Table 1. Summary of the simulations discussed in this paper. The 
number of SPH particles for the simulations of the first section of 
the table is 4 X 10 4 , whereas for the last simulation is 4 X 10 s . 



Run 


Aftot 


(M & ) 


Composition 


Ro (R Q ) 


t (s) 


1 


0.4- 


fl.2 


He/ONe 


0.040 


180 


2 


0.4- 


f0.4 


He/He 


0.042 


600 


3 


0.6- 


f0.6 


CO/CO 


0.041 


180 


4 


0.6- 


fl.O 


CO/CO 


0.038 


725 


5 


0.6- 


f0.8 


CO/CO 


0.033 


600 


6 


0.8- 


fl.O 


CO/CO 


0.028 


1000 


7 


0.6- 


f0.6 


CO/CO 


0.040 


90 



Moreover, the two white dwarfs preserve their initial spher- 
ical symmetry. Hence, and according to Eqs. (II 2> and il'A) . 
we should expect to obtain dimensionless strains which are 
sinusoidal functions with a period equal to half of the or- 
bital period. As it can be seen from Figure |2] the numerical 
solution matches very well the theoretical one. In particular, 
both the amplitude and the frequency (v ~ 0.025 Hz) show 
an excellent agreement between theory and simulations. To 
further illustrate this overall excellent agreement, in Fig- 
ure [|3 we show the residuals between the theoretical solution 
and the gravitational wave emission of the close white dwarf 
binary system in circular orbit. Note that the scale in this 
case is one order of magnitude smaller than that of Figure|5| 
Hence, we conclude that our SPH simulations can accurately 
compute the emission of gravitational waves from coalescing 
white dwarfs. 



5 GRAVITATIONAL WAVE EMISSION 

We have computed the emission of gravitational waves re- 
sulting from the merging of several close white dwarf bi- 
nary systems. In particular, the radiation of gravitational 
waves from a 0.4+0.4 Mq He-He, a 0.4+1.2 Mq He-ONe, 
a 0.6+0.8, a 0.6+1.0 and a 0.8+1.0 Mq CO-CO was com- 
puted. All the simulations performed so far are listed in 
Table 1, where the mass of both components of the binary 
system, their respective composition, the initial separation 
(-Ro) and the simulation time can be found. For the sake of 
conciseness we will only discuss in some detail the results 
of the 0.4+0.4 Mq He-He system and of the 0.8+1.0 Mq 
CO-CO merger. In all the cases the initial separation was 
larger than the corresponding Roche lobe radius of the less 
massive component. For instance, for the case in which a 
0.4+0.4 Mq He-He merger is considered the initial sepa- 
ration was ~ 0.042 Rq, and in the case of a 0.8+1.0 Mq 
CO-CO system the initial separation was ~ 0.027 Rq. In- 
stead of computing self-consistently the chirping phase we 
have chosen to add a very small artificial radial acceleration 
term which decreases the separation of both components un- 
til the last stable orbit is reached. This acceleration term is 
proportional to the velocity, never amounts to more than 
a 5% of the real orbital acceleration and is added once the 
stars have already done a full orbit, then we let the system 
relax during another orbit and the whole procedure is re- 
peated again until the secondary fills its Roche lobe. Once 
the secondary fills its Roche lobe this acceleration term is 
suppressed and the system is allowed to evolve freely. Nev- 
ertheless, we have checked (see §6) that the amplitude of 
the gravitational waves during the initial phase of the co- 
alescence agrees with that of the chirping phase. Finally, 
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Figure 4. Temporal evolution of the 0.8+1.0 Mq CO-CO close white dwarf binary system during the most important phases of the 
merger. The SPH particles have been projected in the orbital plane. See text for additional details. 



it is important to mention here that in all the simulations 
studied in the present work the number of particles for each 
white dwarf is 2 x 10 4 in all cases. However, and in order 
to check the sensitivity of our results to the number of par- 
ticles, we have run an additional simulation in which the 
number of particles was significantly increased to 2 x 10 5 for 
each star, this simulation is listed in the last entry of table 
1 and discussed below. 

In Fig. 2] we show the temporal evolution of the posi- 
tions of the SPH particles projected on to the orbital plane 
as a function of time for the 0.8 + 1.0 Mq system. The ini- 
tial configuration of the two stars was completely spheri- 
cally symmetric. After some time the secondary is tidally 
deformed (top left panel) and begins to overflow its Roche 
lobe (top right panel). As a consequence, an accretion stream 
forms (middle left panel). This accretion stream is directed 



towards the primary, forming an arm. The particles flow- 
ing from the secondary onto the primary are redistributed 
over the surface of the primary (middle right panel) and the 
arm twists as time increases (bottom left panel), leading 
to the formation of a heavy accretion disk with cylindri- 
cal symmetry (bottom right panel). The whole process lasts 
for about 3.8 minutes. The merging process can be under- 
stood in terms of the positive feedback experienced by the 
secondary: as the coalescence proceeds the secondary looses 
mass and, thus, becomes less dense and expands leading 
to an enhanced mass-loss rate which, in turn, leads to a 
decrease of the average density of the secondary. Very few 
particles achieve velocities larger than the escape velocity 
and, hence, a very small fraction of the total mass is ejected 
from the system (about 6.5 x 10 -3 Mq in this case). An im- 
portant feature of the simulation is that the accretion disk 
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Figure 5. Same as figure|l]fbr the 0.6+0.6 Mq CO-CO close white dwarf binary system in which 2 X 10 5 SPH particles were used. Only 
1 out of 4 particles are shown. 



is supported by its own rotational velocity. However, it is 
important to realize that in the final configuration a weak 
spiral pattern can still be found. This pattern should become 
less and less apparent as time increases, reaching cylindri- 
cal symmetry during the very late stages of the simulation. 
However, following the long term evolution of this heavy ac- 
cretion disk would require a heavy computational load which 
is well beyond our current possibilities. Additionally, the ac- 
cretion rate onto the primary becomes negligible during this 
phase. Hence, our final configuration consists in a central ro- 
tating spherically symmetric compact star surrounded by a 
keplerian and almost cylindrically symmetric accretion disk. 
As discussed in Guerrero et al. (2004) the rotational veloc- 
ity of the central star has been problably overestimated due 
to the large shear introduced by SPH methods. We have 
used the artificial viscosity of Balsara (1995), which does 



not produce an excessive shear and, consequently, reduces 
somewhat (but not completely) these problems. Therefore, 
some properties of the merged object could be affected by 
the excess of rotation of the primary. However, since even in 
this case the star preserves spherical symmetry we consider 
that the calculations described below provide a good ap- 
proximation to the emission of gravitational waves. All the 
cases studied present more or less the same features except 
those in which two white dwarfs of equal mass are involved. 
In such a case the final configuration is a single spheroidal 
central object — see, for instance, Fig. 5 of Guerrero et al. 
(2004). 

In figure |S| the coalescence of a binary system of 
0.6+0.6 Mq in which each of the stars was modelled us- 
ing 2 x 10 particles is displayed. As in figure [I] we have 
chosen to represent the temporal evolution of the SPH par- 
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Figure 6. Gravitational wave emission from the merger of a 
0.8+1.0 Mq CO-CO close white dwarf binary system. The di- 
mensionless strains h± and h x are measured in units of 10 — 22 . 
The leftmost thin dotted line corresponds to the time at which 
the two white dwarfs start to merge and the rightmost thin dotted 
line corresponds to the time at which an approximate cylindri- 
cal configuration has been already achieved. Again, the source is 
assumed to be at a distance of 10 kpc. 



tides projected in the orbital plane. Note, however, that in 
this case only one out of four particles has been represented. 
In this case, however, the simulation only covers times larger 
than that at which the last stable orbit of the system occurs. 
As it can be seen, the results are essentially the same and, 
thus, we are confident in the main results and general trends 
of our numerical simulations. We will, however, come back 
later to this issue at the end of this section, when discussing 
the emission of gravitational waves. 

An example of our results is shown in figures |B] and 
where the dimensionless strains h+ and /i x as a function of 
time for different inclinations are respectively shown for the 
0.8+1.0 M Q CO-CO and the 0.4+0.4 M Q He-He systems. 
The beginning and the final time of the merging itself are 
shown in both figures as thin dotted lines. In figure|5]it can 
be seen that before the coalescence proceeds the emission of 
gravitational waves still has a sinusoidal pattern, but with 
an increasing frequency. That is, the close white dwarf bi- 
nary system chirps as a consequence of the spiral trajectory 
of the stars towards the center of mass. Note that for 
i — 7r/2 is zero because the orbital plane is parallel to the 
line of sight. When the two white dwarfs start to coalesce, 
the amplitude of the dimensionless strains somehow increase 
first, but only during the first (and most violent) stage of the 
merger. This corresponds to the phase in which a spiral arm 
is formed. After the second maximum is achieved the am- 
plitude decreases dramatically. In fact, only two more clear 
maxima can be apparently distinguished before the system 
reaches its final configuration. It is interesting to note that 
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Figure 7. Same as Fig.|S]but for the case of a 0.4+0.4 Mq He-He 
close white dwarf binary system. 

once the merger has already finished one of the dimension- 
less strains still is significant, h+ ~ —5 - 10 -22 at d = 10 kpc 
for i = 7t/2. This residual emission is due to the inhomegeni- 
eties of the accretion disk previously discussed. Each one of 
the very small maxima appearing at very late times corre- 
sponds to successive crossings of the edge of the weak spiral 
arm in front of the line of sight. Note as well that this resid- 
ual emission tends to disappear assymptotically, as a con- 
sequence of the ongoing rehomogeneization of the accretion 
disk. It is important to realize that these inhomogeneities 
could be either an artifact due to the resolution used in our 
SPH simulations or a consequence of the adopted artificial 
viscosity, since it is well known that the artificial viscosity 
of Balsara (1995) induces a considerable shear viscosity. 

In figure Q it can be seen that although during the first 
part of simulation the same chirping pattern is found, once 
the merger proceeds the gravitational wave signal suddenly 
disappears on a short time scale, comparable to the orbital 
period. This is due to the fact that in this case the two stars 
have equal masses and, therefore, there is not a prominent 
spiral accretion stream. Instead the two components of the 
binary system are disrupted around the center of masses, 
where a shocked region forms as a consequence of the im- 
pact between two streams. The typical Mach numbers in 
the shocked region are Ma ~ 1. The streams are slightly 
asymmetric, depending on the orbital phase at which they 
form. By the end of the simulation — see again Fig. 5 of 
Guerrero et al. (2004) — both streams become entangled 
and the final configuration of the resulting object consists 
in a central shocked region surrounded by a less dense ro- 
tating spheroid, in which a certain degree of asymmetry is 
still present. Hence, the coalescing process does not present 
a rather symmetric behavior during the initial phases and, 
consequently, we see three consecutive maxima. After this, 
once the streams become entangled the emission of grav- 
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Figure 8. Modulus of the Fast Fourier Transform of the adimen- 
sional strains h+ and h x of Fig. 



Figure 9. Modulus of the Fast Fourier Transform of the adimen- 
sional strains h+ and h x of Fig. 171 



itational waves is heavily suppressed and, by the end of 
the merging process, the emission of gravitational waves is 
negligible. However, note as well that in this case a small 
residual emission is also observed in the dimensionless strain 
h+ = —0.8 ■ 10 -22 at d — 10 kpc for i = tv/2, although con- 
siderably smaller than in the previous case. This is again due 
to our limited computational resources. The impossibility of 
following the very late phases of the coalescence episode does 
not allow us to compute the rehomogeneization of the exter- 
nal spheroid and, hence, to compute the long-term behavior 
of h+ accurately. 

The same information is displayed in figures |H] and |3] 
but in a different format. As can be seen there, the dominant 
frequency is given by the orbital period. However, as the 
components of the binary system approach each other, the 
dominant frequency is shifted to larger values and, during 
the merger, high frequencies show up, although with very 
low amplitudes. 

Finally, and in order to check the sensitivity of our 
results to the resolution of the SPH simulations described 
above we have computed the emission of gravitational waves 
for both our run "7", in which 2 x 10 particles were used 
for each star, and our run number 3, in which a resolution 
10 times poorer was adopted. The results are displayed in 
Fig. E3 As in the previous figures the vertical thin lines de- 
note the moment at which the last stable orbit is achieved. 
Note, nevertheless, that the high resolution simulation was 
started, as previously mentioned, when the system was al- 
ready at the last stable orbit. That is the reason why the 
chirping phase does not appear. Consequently, the time ori- 
gin of this simulation has been shifted accordingly to match 
that of the low resolution simulation. As it can be seen the 
emission of gravitational waves is very similar in both cases, 
and, hence, our results are robust. 



6 DETECTABILITY 

In order to check whether or not LISA would be able to de- 
tect a close white dwarf binary system we have proceeded 
as follows. We have already shown that the most promi- 
nent feature of the emitted signal is its sudden disappear- 
ance in a couple of orbital periods and that the gravita- 
tional wave emission during the coalescence phase does not 
increase noticeably. Hence, the gravitational wave emission 
is dominated by the chirping phase. Hence, we have assumed 
that the orbital separation of the two white dwarfs is exactly 
that of our binary system when mass transfer starts. We 
have done so because, as explained before, we have added a 
small artificial acceleration term to the initial configuration 
in order to avoid an excessive computational demand at the 
very beginning of our simulations. This acceleration term is 
suppressed once the secondary begins to transfer mass onto 
the primary. Note, however, that the mass transfer starts 
when the secondary fills its Roche lobe and, consequently, 
this orbital separation is physically sound. We have further 
assumed that the integration time of LISA will be one year. 
We have checked that during this period the variation of the 
orbital separation is negligible (see also figures |S| and 0. Of 
course, should the integration time be smaller the signal-to- 
noise ratio derived below would smaller. Hence, our results 
should be regarded as an upper limit. 

The signal-to-noise ratio, rj, is given by 



2 

V = 



h 2 (u>) duj 
S(u>) 2^ 



(14) 



where S(ui) = Sh(w)r is the sensibility of LISA, r is the 
integration period, and h(u>) is the Fourier Transform of 
the dimensionless strain. It can be easily shown that for 
a monocromatic gravitational wave n = h{u) / S]/ 2 (uj). The 
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Table 2. Maximum distance at which LISA will detect the close white dwarf binary systems discussed in this paper, see text for details. 
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Figure 10. A comparison of the computed emission of gravi- 
tational waves when the resolution of our SPH simulations is 
changed. The left panels show the emission of gravitational waves 
for a simulation in which a low resolution (2 X 10 4 SPH particles) 
was used. The right panels show the same quantities when a high 
resolution was adopted (2 X 10 5 SPH particles). See text for fur- 
ther details. 



maximum distance, d ma x, at which LISA would be able to 
detect a close white dwarf binary system is then: 



dmax ~ 17 ( - 



1\, 



M \5/3 

MqJ ll miiz ' 



n2/3/ 10-^ 



kpc (15) 



where vq is the frequency of gravitational wave, and M is 
the chirping mass. 

(rf) 3/5 (16) 



M 



being jj, — mimz/Mtat the reduced mass, and M to t = mi + 
m,2 the total mass. 

In order to evaluate the maximum distance at which 
LISA would be able to detect the close white dwarf binary 
systems studied here we have adopted r\ = 5. We have fur- 
thermore used the integrated sensibility of LISA as obtained 
from http: //www. srl . caltech.edu/~shane/ sensitivity 



Figure 11. A comparison of the signal produced by the close 
white dwarf binary systems studied here, when a distance of 
10 kpc is adopted, with the spectral distribution of noise of LISA 
for a one year integration period. 



The results are given in Table 2, where the frequency of 
the close white dwarf binary system when the secondary 
overflows its Roche lobe, the chirping mass, the maximum 
distance at which LISA would detect them for a signal-to- 
noise ratio of 5 and one year integration, and the signal-to- 
noise ratio at 10 kpc for one year integration are shown. Also 
shown in Table 2 are the peak amplitude at 10 kpc of the di- 
mensionless strain (in units of 10 -22 ) and the total radiated 
energy in the form of gravitational waves. We stress that 
if the merger occurs during the one year integration period 
the signal-to-noise ratio would be smaller. It is worth notic- 
ing as well that the energy radiated away from the binary 
system during the coalescence in the form of gravitational 
waves is of the order of ~ 10 41 erg, much smaller than the 
total energy of the system and, hence, totally negligible in 
the energy budget, so no back reaction should be taken into 
account and the procedure used here is robust. Note as well, 
that for the case of 0.6 + 0.6 Mq system, in which two res- 
olutions were used the maximum distance at which LISA 
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would be able to detect the coalescence is very similar for 
both the high resolution simulation, 12.8 pc, and the low 
resolution simulation, 13.0 pc. The energy released during 
the coalescence is very similar as well in both cases, 0.63 and 
0.64 x 10 41 erg, respectively. 

Finally, in Fig. El we compare the signal produced by 
the close white dwarf binary systems studied in this paper, 
when a distance of 10 kpc is adopted, with the spectral dis- 
tribution of noise of LISA for a one year integration period. 
As it can be seen, all of them will be eventually detected, 
at different signal-to-noise ratios, ranging from ~ 5.0 for the 
0.4 + 0.4 Mq system to ~ 16.6 for the 0.8 + 1.0 Mq system. 



7 DISCUSSION AND CONCLUSIONS 

We have computed the emission of gravitational waves of 
merging white dwarf binaries, for a wide range of masses 
and compositions of the components of the binary system. 
For that purpose we have used a SPH code which allowed 
us to follow the temporal evolution of the coalescing white 
dwarfs. We have shown that the most noticeable feature of 
the emitted signal is a sudden disappearance of the gravita- 
tional strains. By contrast the chirping phase will be easily 
detectable by future space-borne interferometers like LISA. 
In fact, it can be said that the most relevant signature of the 
merger will be the absence of any signature and the sudden 
disappearance of the source. Since the frequency coverage 
of LISA will range from 10 _1 to 10 -4 Hz, the detection 
of chirping close white dwarf binary systems is guaranteed 
(Farmer & Phinney 2003; Nelemans et al. 2001a). Moreover, 
at frequencies of ~ 3 • 10~ 3 Hz most close white dwarf bi- 
nary systems will be spectrally resolved (Cornish & Larson 
2003). Typically, LISA will be able to detect ~ 3000 bi- 
naries at these frequencies (Seto 2002). Additonally, it has 
been recently shown (Cooray, Farmer & Seto 2004) that the 
typical error box of LISA for these kind of systems will be 
Sfl ~ 5 deg 2 , and that, given that many of the detected 
sources will be eclipsing binaries with a period equal to that 
of the gravitational waves, optical follow-up campaigns will 
allow us to further constrain the location of the sources. 
Hence, by combining optical observations and gravitational 
wave data, and taking into account that the gravitational 
wave signal for such mergers suddenly vanishes — or, equiva- 
lently, that during the integration period the signal-to-noise 
ratio stops growing — we should be able to gain insight 
into the physics of merging and to obtain precise informa- 
tion of the properties of the progenitor systems. Moreover, 
it should be taken into account that for a typical merger 
— namely, the 0.6 + 0.8 Mq case — the volume accesible 
to LISA is Vlisa ~ 1.2 x 10 14 pc 3 . Since the volume of the 
Galaxy is V ~ 3 x 10 11 pc 3 , LISA would be able to detect 
all the mergers occuring in our Galaxy during its operation 
period. However, the typical rate of white dwarf mergers is 
r ~ 8.3 x 10~ 3 yr _1 (Nelemans 2003) and, hence, although 
there is an uncertainty of a factor of 5 in the rate of white 
dwarf mergers, the expected detection rate is consequently 
small. 
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